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Abstract. We review some of the recent developments in two dimensional statis- 
tical mechanics in which corner transfer matrices provide the vital link between 
the physical system and the representation theory of quantum affine algebras. 
This opens many new possibilities, because the eigenstates may be described us- 
ing the properties of g-vertex operators. 
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1. Introduction 

Transfer matrices have been a central object for the study of two-dimensional 
statistical mechanics systems for more than a half century. Originally the row transfer 
matrix was defined: it encapsulates the statistical information (via the Boltzmann 
weights) pertaining to the contribution from a single row of a regular lattice. Row 
transfer matrices are extremely complicated; for an Ising or six-vertex model on a 
square lattice of width L the row transfer matrix has dimension 2^. Onsager's solution 
of the Ising model [1] used representation theory of finite- dimensional Lie algebras to 
transform the generators of the row transfer matrix to a tractable form. The next 
several major advances — for example the solution of the one-dimensional Bose gas 
[2] , the six [3] and eight [4] vertex models — brought a totally different approach to 
the forefront, the Bethe Ansatz. But developments in the last decade have restored 
the role of representation theory. The quantum inverse scattering method [5,6], an 
algebraic implementation of the Bethe Ansatz and Yang-Baxter equations, gave birth 
to the theory of quantum groups [7,8]. Representation theory of infinite dimensional 
Lie algebras became paramount in conformal field theory [9,10]. 

An independent development was the invention of the corner transfer matrix 
(CTM) [11,12]. It is the partition function for a whole corner of a regular lattice and 
might therefore be expected to be even less tractable than the row transfer matrix. 
Surprisingly, this is not so: CTMs have an extremely simple eigenvalue spectrum (in- 
teger powers of a single parameter). CTMs proved to be an effective method for the 
evaluation of one-point functions [12-14]. In [14] configuration sums appeared which 
were identified as characters of Virasoro algebras [15,16]. These mysterious connec- 
tions led to the discovery of hierarchies of solvable lattice models in two dimensions 
[17] and of many beautiful connections with infinite dimensional Lie algebras [18-20]. 
Furthermore, the multiplicities of the CTM eigenvalues are equal to the weight-space 
multiplicities of irreducible highest weight representations of certain affine Lie alge- 
bras [20] , and the actual configurations used to label the CTM eigenvectors also label 
the crystal base vectors of the corresponding representations [21-24]. These were vital 
clues for finding the links with representation theory. But the crucial point was to 
identify the appropriate algebras and their representations, and the role of the row 
transfer matrices and CTMs in these algebras [25,26]. 

In one of the original papers on CTMs [11], Baxter observed that there is no 
Ising-like reduction in the six and eight vertex models. He wrote 

A rather ambitious hope is that by examining the CTMs we may stumble on 

such a group, that the solution of the models may thereby be simplified ... 

In [26] was given a new scheme for solving the six-vertex model, in the antiferromag- 
netic regime, using the newly discovered quantum affine symmetry of the CTM. The 
approach of that paper has been extended to higher spin chains [27], to the higher 
rank case [28] and to the ABF models of Andrews, Baxter and Forrester [29]. The 
scheme has already been used in some of those models to give general expressions for 
the 71-point correlation functions [30]: it also provides a means to derive g-difference 
equations for them [30,31]. The promise of this new approach is self-evident. Many 
hard-won results of Bethe Ansatz calculations have been recovered, important new re- 
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suits have been obtained, and there is much research in progress, since the calculation 
of correlation functions is a long-standing problem. 

The six- vertex model is related to the quantum affine algebra Uq(sl2) of Drinfeld 
and Jimbo [8,32], since its Boltzmann weights are the simplest example of an R- 
matrix. It is a paradigm for the application of the quantum affine symmetry to the 
representation of CTMs. Briefly, the magic of the six-vertex CTM arises from the 
property that its eigenvectors provide level- 1 highest weight modules of the quantum 
affine algebra Uq(sl2), while the CTM acts as a derivation operator in the algebra. 
Similar statements may be made in the other cases (higher spin, higher rank). Stated 
more simply, the Chevalley generators of the quantum affine algebra are raising and 
lowering operators when acting on the eigenvectors of the CTM, and the importance of 
this cannot be over-stated. It is precisely the kind of property which makes the Ising 
model so tractable, although the fermion algebra is much simpler than the quantum 
affine algebra with which we deal here. 

The purpose of this article is to give an introductory account which is as ele- 
mentary as possible. It is based on talks which have been given by the author. For 
simplicity we consider in detail only the six-vertex model on which the original ideas 
were developed. In addition, the focus is on the role of CTMs; consequently we omit 
any discussion of the construction and properties of creation/annihilation operators 
for the XXZ Hamiltonian, although this was a major goal in [26] . There is a growing 
literature concerned with the extension of the methods to various statistical mechanics 
systems and with the use of g-difference equations which are an essential ingredient in 
the theory of g-vertex operators. The six-vertex model uses the level-1 modules of the 
rank-1 quantum affine algebra Uq(sl2)- Higher-spin and higher-rank cases are consid- 
ered in [27,28,33]. The Boltzmann weights of all these models form the i?-matrices 
which intertwine tensor products of finite-dimensional representations. The relation 
is more complicated for the ABF models: the Boltzmann weights are the connection 
coefficients which intertwine vertex operators, themselves intertwiners [34]. The pro- 
cedure for the ABF models [29] is somewhat analogous to the Goddard-Kent-Olive 
coset construction of the discrete minimal series of Virasoro algebras [35]. For the 
many details of all these developments the reader should consult the cited literature. 

2. Transfer matrices 

2.1 Definition of the model. The six-vertex model is defined on a regular lattice, each 
vertex of which is at the intersection of two edges, which carry the physical spin 
variables ±. The energy Ec of a configuration of the system derives is the sum of 
the energy assigned to each vertex. From this additivity, the Boltzmann weight 
exp{—E(./kT) of a configuration is the product exp(— £'v//cT) of vertex weights 
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The six- vertex model is "Z-invariant" in the sense of Baxter [36]. For this reason we 
attach directions and rapidity variables Ua-y, ups to the lines of the lattice, so that 
W{a,P,^, 5; u) has as its rapidity their difference. The convention is that negating a 
direction negates the corresponding Ua^y or ups- 

Only six of the Boltzmann weights are non-zero. In the notation of [37] they are 
a = psinh(A — ti), b = psinhu, c = psinhA, where p is a normalisation factor and 
< « < A in the physical regime. In anticipation of the connection with the quantum 
affine algebra, we introduce new variables by 

q = - exp(-A), C = exp(w), A = ^'^^ . 
Also, we normalise the weights by the partition function per site, k{(): 
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2.2 Properties of vertex weights. The most important property of the weights is that 
they satisfy the Yang-Baxter equations, which have the well-known pictorial repre- 
sentation 




This is the origin of Z-invariance. Summation is implied over all internal lines and 
this sum is invariant under rearrangement of the lines, provided only that all vertices 
have the same value of A and the boundary arrangement is fixed. Because of the 
locality of the Yang-Baxter equations, this holds for arbitrarily large systems. 
Two other important properties are unitarity 
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and crossing symmetry 
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The partition function per site may be obtained from these by the standard 
inversion trick [37]. We sketch the argument. Unitarity and crossing symmetry imply 
the relations 

«(C) <C') = 1, n{~q-'C') = ^(TT^ '^(0- 

This gives a ^-difference equation and the unique solution for a function analytic in 
the annulus < |C| < q~^, and normalised to = 1, is 

"^^^ AC-'q';q')oo{eq';q')oo' 



where 



n=0 

2.3 Corner transfer matrix. CTMs are defined as the partition function of a quadrant 
of the lattice. The configurations of the edges of the quadrant label the rows and 
columns of the CTM whilst the outside boundary, which is part of the boundary of 
the total system, is fixed so as to select the required ground state. There is a preferred 
diagonal direction (SE-NW), which is along the boundary of A{u), and two ground 
states labelled z = 0, 1, related by a lattice shift of one site. 

CTMs are defined in the infinite lattice limit in a normalised form [37]. A small 
part of the construction of two CTMs, A{u), B{v) looks like: 
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We consistently use the notation for matrix elements of an operator L : V ^ W that 
Lva = L^wp: the upper index labels the row. 
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Recall Baxter's arguments from which the properties of the CTMs stem. He 
considers the product of two CTMs A{u), B{y). The result is the partition function 
of an infinite half plane and its value is a matrix element of an arbitrarily high power 
of the row transfer matrix. This depends only on the maximal eigenvector, which 
in turn is a function only of the rapidity difference u — v (Yang-Baxter equations). 
This implies B{y)A{u) = X{u — v). The group property A{v)A{v) = A{u + v) is the 
most important conclusion to be drawn. It follows from the crossing symmetry of the 
Boltzmann weights and spin reversal. The result is 

where 

—2 °° 

* n=l 

and R is the operation of spin reversal at every site. 



3. The quantum affine algebra 

3.1 Notation. In this section we define our notations for Uq[si2) and review those 
properties needed for our discussion of the six- vertex model. t/g(sl2) is generated by 
Ci, /i, ti = {i — 0, 1), which satisfy the defining relations [38] 

tiCj = Q Cjti, tifj ~ Q fjtii 

titj — tjti, [e^, fj^ — ^^'^'q g— 1 ' 

efej - [3]e^eje, + [3]e,eje2 - Cjcf = 0, 

fffj - [3]/f + [3]/./,/f - fjf! = 0, ^ 

Here [n] = (g" — q~^)/{q — q~^), and Aij is the generalised Cartan matrix for the 
affine Lie algebra 5I2 

'2 -2' 
-2 2 



Aij — 



To obtain the full quantum affine algebra Uq(sl2) one must add the generator g'^. d 
is the derivation, and it satisfies 

[d, Si] = Si^oCi, [d, fi] = -Sififi. (1) 

The co-multiplication for the Hopf algebra structure is defined by 

A(ei) = ei®id + ti® e^, 
^{fi) = fi®tr^ +id®fi, 
A{ti) = ti®ti. 
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and the formula for the antipode is 

a{ei) = -t~^ei, a{fi) = -fiU, a{ti) = t~^. 

3.2 The R-matrix. The simplest spin half Uq(^si2) module V over the field Q{q) has 
basis vectors v+ = vq, V- = vi, and is equipped with the actions 

eiV+=0, fiv+^v-, tiv+ = qv+, 

eiv- = f+, fiv- = 0, tiv- = q~^v-. 

From this one may construct a Uq(^sl2) module by introducing the formal variable z. 
The basis vectors of this infinite-dimensional module are {z"^ ^ v±} but we simply 
write {z'^v±}. That is, one equips = Q,[z, z~^] (8) V with the Uq(^sl2) module action 

a I— > z—. 

eii-^ei, /I'-^/i, ht-^ti, dz 

Vz is also a Uqisi-^) module: for this one simply regards 2; as a fixed non-zero complex 
number and drops the grading operator d (evaluation representation). 

Define the operator R{z) by 

R{z)v± ®v± = v±®v±^ 



R(z)v+ (g) f_ = 7^qv+ ®v- + T^zv- (g) v+, 

1 — q^z 1 — q^z 



1 — q^ 1 — z 

R(z)v- = 5- v+ ®v- + T^qv- ® v+. 

1 — q^z 1 — q^z 

Further, let P be the transposition operator: P-{vi<S)V2) = {v2^vi). Then R{zi / Z2) = 
PR{zi/z2) : Vzj^ ®Vz2 Vz2®Vz-^ is an intertwiner of (3(2) modules. Its entries also 
provide the Boltzmann weights for the six- vertex model with a different normalisation 
and a change of gauge. Writing z = 

W{i,3,k,l;0 = K{z)-^g-^^R^]{z), 

with ^fjfc = 1, i = k, and g±^i = C^^. The gauge factor is the well-known term which 
makes finite spin chains [/g(s[2)-invariant. It simply adds the total spin to Kxxz 
and changes its role as a grading operator. Since it does not affect the physics, we 
henceforth use the Boltzmann weights 

W{i,j,k,l;z) = gikW{i,j,k,l;C)- 
This is consistent with the cited works [26-29,33]. 

3.3 Dual modules. The dual module construction is intimately connected with the 
crossing symmetry and is required for representing local operators. It goes as follows. 
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Given a linear space V its dual V* is the set of linear maps V* : V ^ Chy the action 
u* 

V I — > {u*,v). Choose the dual basis by the canonical pairing {v^Vk) = Sjk- V* 
becomes a dual module under the action 

{XU* ^v) = {u* ^(j){x)v)^ e t/q(sl2), 

where (/> is any anti-homomorphism of the algebra {(plxy) — 4>{y)4>{x)). Using the 
antipode for 0, the action of Uq(^sl2) on V* is 

eiv^ = -q v_, fiv^ = 0, tiv^ = q v^, 

eiv*_ = 0, /iv* = -qv"i, tiv*_ = qvt, 

eo ^ {zq~^)fi, fo ^ {zq~'^)~^ei, to ^ t^^. 

This gives the isomorphism 

C : V* — > Vgq-2, ^ V-, I— > —qv+. 

The operator R{z) was defined to act on T^i '^Vz2 so that PR{z) is an intertwiner. 
We also require the equivalent operator R*{z) acting on V*^ ® Vg^. It is given as 
R*{z) = {C ^id)-^R{zq-^){C ^id). Explicitly, 

n*/ \ * q — q~^z ^ — * 

R {z)vX (g) v+ = - — (8) v+ + - — - — v*_®v-, 

1 — z 1 — z 

R {z)v_ (g) f _ = zv, (g) u+ H v_ ^ 

1 — z 1 — z 



3.4 Notation. Our notations for sl2 are as follows. The Cartan subalgebra is 
spanned by {ho^hi^d} and ao, ai are the roots. They are related to the funda- 
mental weights by ao = 2Aq — 2Ai + d, ai = 2Ai — 2Ao; we also write p = Aq -|- Ai. 
The invariant form on ^* is given by 2(Ai,Aj) = ^i^Jij, (Ai,5) = 1, (5, 5) = 0. The 
weight lattice and its dual are P = ZAq © ZAi ® Z6 and P* = Zho ® Zhi ® Zd, with 
(Aj, hj) = Sij, {Ai, d) = 0, {5, hi) = 0, {6,d) = 1. We identify P* with a subset of P 
via ( , ), so that ctj = hi and 

2p = 4d + /ii. (2) 

This plays an important role in relating Kxxz and i^cTM with p and d. 

There are two irreducible infinite dimensional level 1 highest weight Uq {sl2) mod- 
ules, V{Aq), V{Ai), corresponding to the fundamental weights Aq, Ai. Their relation 
to the eigenvectors of the CTM will be discussed in the next section. 
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4. The CTM and its eigenvectors 

In [25] Foda and Miwa identified a six-vertex CTM generator witli ttie derivation 
d of Uq (SI2) 7 and its eigenvectors witli tlie weiglit vectors of the level- 1 modules. There 
are two bare ground states, {■■■,—,+, — , -|-) and {■■■,+,—, +, —). They defined the 
action of the algebra on the sets of paths Vi, (z = 0, 1), which differ from these states 
at only a finite number of places. The problem is to show that everything works 
consistently in each order of perturbation. Convincing evidence was given for this, 
through extensive analytic computations to quite high powers of q, but the main 
result remained a conjecture. Moreover the g-expansions about these bare states are 
divergent, and an infinite renormalisation is involved. 

In the next two subsections we explain briefly the approach of Davies [39] , which 
uses finite size truncation together with the Kashiwara's crystal base theory [21,22] to 
prove the necesary identifications. This method generalises quite easily to the higher 
level and higher rank cases. 

4.1 Identification of derivation. Kxxz is defined using the Boltzmann weights W, but 
it is convenient to introduce another generator Kctm defined from W. It is built from 
two-site operators Hj-. 

obtained by expansion of R about z = 1. The subscript j indicates that Hj operates 
at the positions j + 1, j oi the chain. For this reason, we attach a second subscript 
to the generators e^, fi, U, writing Cij, fij, tij. The truncated CTM generator is 

N-l 

-^CTM ~ 2-^ J^J- 

It acts on the A'"-fold tensor product F®^. Write e^, fi, U for the A'"-fold iterated 
co-product of the generators, 

N 

Ci = ^ ti^N ® • • • (8) U^j+i ® Bij (8) idj-i (8) • • • (8) idi, 

N 

= 5^ idjv <8) • • • <8) idj+i (g) fij (g) t~j_;^ (g) • • • <8) t~l , 



'i,N (X) • • • 



We need their commutation relations with K^-^^. For ei, /i this is trivial, 
R{x,y) intertwines the action of the subalgebra Uq(^sl2) which does not depend on 

the spectral parameter. If follows immediately that -P^ctm Uq(5l2) invariant: 

[4™,ei]=0, [Kg:)^,A] = 0. 
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The intertwining condition for cq, fo gives the commutation relations [39] 
[-^CTM' ^o] = eo - A^eo,Ar (g) idA^-i ■■■(g) idi, 

This is the vital result. In the infinite N limit the boundary terms may be neglected, 

so that KcTM = liniiv^oo(-^cTM ~ ^n) acts as the derivation in the algebra Uq(^3i2), 
verifying the relations (1). i?iv is a renormalisation constant and plays no further 
role. 

4.2 Identification of modules. Kxxz and i^cTM act formally on a semi-infinite tensor 
product of spin-half modules, as does the quantum affine algebra. It must be em- 
phasised that this limit is only formally defined. When the eigenstates are developed 
either as ^2-sequences or as g-expansions, one obtains divergent expressions. The sit- 
uation is similar to the Ising model CTMs [40,41]. Having truncated the chain at a 
finite length, the action of the algebra is well-defined. K(yYM renormalised by its 
ground state eigenvalue. It is shown in [39] that the eigenvectors with eigenvalues 
up to M carry the ?7g(sl2) action to within error terms of order qi2-^v-2M-i_ ^j^g 
boundary terms carry factors of the order Nq'^^^ , which is why they may be ne- 
glected in the infinite limit. It is also shown there, using the theory of the crystal 
base, that the eigenvectors are in one-to-one correspondence with the weight vectors 
of the modules V{Ai). 

This method of finite size approximation had already been noted for the crystal 
base [24], so it is natural to extend it to g 7^ 0. The meaning of the approximation 
is quite precise, and illustrates the non-uniform nature of the convergence. For any 
fixed integer M, the eigenstates corresponding to the grading levels < M approximate 
weight vectors up to that grading level to within errors of order ]Vg^-^~'^^~^ as 
N 00. But if one selects some order of approximation, and then increases M and N 
together, the approximated eigenvectors become an increasingly small fraction of the 
totality of eigenvectors of the truncated CTM, since 2^/2^^ increases exponentially 
fast. What matters is that g-expansions of physical quantities obtained from the CTM 
and from representation theory will agree to every order in q. 

Without further ado, we shall identify the eigenvectors of the CTM with the 
weight vectors of the highest weight modules, and the generator Kqtu with the 
derivation operator d. Examination of the gauge term in W shows that i^xxz = 
2-f^cTM + S., 2S — cr|. Using the connections (2) leads to the identification 

Kxxz = 2d + hi/2 = p. 



4.3 Local properties. Let L be a local operator, acting at the sites 1, ■••,n. The 
expectation value (L) is calculated by the standard CTM method [37] as a quotient 
of traces: 

_ trjLABCD) _ irv^K^)W) , . 

^ ^ ir{ABCD) trviK,){q^f>) " ^ ' 
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Here A, B, C and D are CTMs for the four quadrants of the system. The difficulty 
is that L must be expressed in a basis of eigenvectors of i^xxz- Without an exphcit 
construction for these, the calculation is restricted to an operator which commutes 
with Kxxz, i-e., a one-point function. For the six- vertex model this L = (— 1)"^ 
and the one-point function is the magnetisation in the "interaction round a face" 
representation [37]. The traces are specialisations of the character formulae: 



Xt{x,y) = trv.(A,)(a; y = y 



Substituting x = q^, y = q into these formulae gives the magnetisation as 



m 



n 



5. Local description of CTM eigenvectors 

5.1 Change of basis. The difficulty of defining the CTM eigenvectors as normalisable 
£2 sequences has already been noted. Representation theory provides a way around 
the problem. Consider the tensor product y(Ai_j) ® V. Formally, we identify this 
with {^'^^2^n) ® Vi, with boundary condition for sector z, but counted from the first 
site. (There is some abuse of notation: Vn here denotes the spin-half module for site 
n.) The product may also be written ((^^t^iKi) and then equated to weight vectors 
of V{Ai). The map between the two pictures is a q-- vertex operator (type I Vertex 
Operators or VO [34]) of Frenkel and Reshetikhin [34]. 



In this diagram crosses represents sites in the chain (copies of V), not vertices in 
transfer matrices. 

VOs are simply homomorphisms — they preserve the algebra structure. They 
may also be interpreted as expansions of the CTM eigenvectors for boundary sector 
i in terms of the state variable at site 1 and the eigenvectors for the boundary sector 
1 — i. This process may be iterated to obtain information about any finite number of 
state variables at sites 1, . . . , n. 

VOs give precisely the information required for the calculation of local properties 
in (3). Moreover, the expansions have convergent domains in q and z. But the most 
crucial property is that the intertwining property of the VOs uniquely determines 
them to within normalisation. Once one has made the identifications of section 4, 
no further physical arguments are required. The introductory comments above, and 
those of section 5.3 below, are quite superfluous in this sense. They are included to 
provide some intuitive understanding of the physical meaning of the theory. 

5.2 Vertex operators. The basic definition of VOs is that they are Ulj(sl2) homomor- 
phisms 

$i;-^^:^(AO^F(Ai-i)®^, 
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We discuss It gives an expansion of weight vectors u e as tensor 

products: 

U=^ Wn,j ® Vj, Wn,j E V{Ai-i). 

3=± 

The intertwining property 

imphes that wt('u) = wt{wn,j) + wt{vj) + nd. Plainly speaking, the expansion is over 
states of all energy whose spins match. The normalisation is fixed by setting the 
coefficient of the term in the highest weight vectors u\. to unity: 

^A- "'^«A, = («Ai_, «) Vi_i) + ■■■, 

VOs may also be defined to preserve the grading operator d: they are the Uq[sl2) 
homomorphisms 

^t-'''{z)■.V{Ai)^V{A^_i)®V,, 

The two forms and ^{z)) are trivially related by inserting or deleting the appro- 
priate powers of z, 

■U = ^ Wn,j <S> Z^Vj. 

3 = ± 

Since the derivation detects the grading before and after the action of the VO, the 
powers of z may be inserted by conjugation with the operator z^: 

From this follows the "homogeneity property" 



Some elementary calculations give the following results: 

^^l^{z)uAo = MAi ^v-- qifiUA,) (^v+ + 0{z), 
^^fiz)uA,=UA,0v+ + Oiz), (4) 

^if{z){fiUA,) = UA, ®V^ + 0{Z). 

They will be used below to fix normalisations. 
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5.3 Physical description. Let Pg'^' be a non normalised semi-infinite shift operator 
with entries 

• • •'^/33,72^/32,7l^/3lJ• 



p■•■,72,7l ,3 

...,/33,/32,/3i 



Since P preserves the U^[sl2) action one may identify an appropriate constant mul- 
tiple with the VO $a'" • Further, we use the relation of d with the CTM to write 



— Ai-i{z)PAi{z~ ). We may even identify the VO as a normalised row 
transfer matrix. This is a generalisation of the boost property for the full transfer 
matrix [42,43]. The product PAi is the partition function of a quadrant with the 
boundary type i along its horizontal edge and 1 — i along its vertical one. It is not a 
CTM, but may be re- written as a different product: 



Ai- 



A^ 



/33- 


73 


f32- 


72 


Pl- 


71 



72 

71 



_ Ps P2 Pi 



P . 

"3 "2 ai 3 



CKs a2 Oil 



The semi-infinite row of vertices corresponds to the vertex operator '^(-2), with 
a normalisation factor 4>i{z). 



This physical interpretation of VOs follows [31]. 

The product of row transfer matrices is commuted using the Yang Baxter equa- 
tions together with the assumption that boundary effects disappear for a system with 
fixed boundary conditions: 



Z2 



Z2 



Zl 



We use this to fix normalisation factors. Write 



(j)l-i{zi)(j)i{z2) 



and take the vacuum expectation. From the definitions, the terms in '^(^i)uAi 
involve only non- negative powers of zi. In order to arrive once more at the highest 
weight vector after the action of ^\^^_.{z2), these terms must pair with a matching 
negative power of Z2- Denote the vacuum expectation as 



*.(^l/^2) = («A, I ^i;r,(^l)*A:-^'' {Z2) I «A,) 
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It involves only non-negative powers of Zi/z2, and will be analytic in Zi/z2 in some 
domain containing the origin. As an element of ^ ® it must be an eigenvector of 
R{zi/ Z2) because of the commutation relation. A simple calculation gives: 

2 

R(z)(v+ ^v- - qv- ® v+) = 5- (f+ ® f_ - qv- ® v+), 

1 — q'^z 

R(z)(v- ®Vj^ — qz ^Vj^ ® v_) = 5 — (v_ — qz ^Vj^ ® 

1 — q'^z 

To leading order this corresponds exactly with (4). Inserting the formula for k{z) and 
using analyticity, 

*o(^) = 14^^ (^+ ® ^- - qv- ® v+\ (5) 

\q z^ q 



00 

— J- — {v-(^v+-qz v+(^v-). 



00 



and 

0i(-s) = 7i-2^^-*"^% Ao = 0, Ai = l/4, 

for some constants 7^. The "conformal factors" Aj, which were included in the defi- 
nition of the VOs in [26], are a shift in the zero point; Ai{z) ~ z'^'^^'- . Such shifts are 
calculated for the Ising model in [41]. 

5.4 Mathematical construction. We have discussed the physical construction of the 
VOs in some detail. However, for the purpose of calculation, the mathematical con- 
struction is preferable. The scale factors (5) are determined unambiguously in [26,28] 
by solving the q-KZ (Knihnik-Zamolodchikov) equation [34] . From this, the commu- 
tation relation for the VOs follows immediately, viz: 



{zq^;q'^)oo{z~^q^]q^)c 

r{z) = 



1 00 



{z-^q^;q^)oc{zq'';q^)oc' 

This confirms our identification of the conformal factors, and the role of the partition 
function, in interpreting the VOs as transfer matrices. 

The charge conjugation isomorphism C relates ^a-"*^ l-^) ^A-"'^!-^) 

and from this the commutation relation emerges in a form needed for the calculation 
of the spontaneous polarisation. It is 

PR%z)^t':*iz^)^^^-^''{z2) = -q-'z'Mzq-')^^^^^^^^ 
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To complete the determination of local operators, we must construct the inverse 
of $^^"'^(^). By definition, it is ^'a1_ normalisation. The construction 

goes as follows. The VOs are algebra homomorphisms, so one may use the canonical 
identification that Hom(L (g) N,M) is isomorphic to Hom(L, M (8) AT*), to equate 

^A^-ivi'^) with (z). Let us expand on this. It means that 

where the pairing is in V* V. One may check that this is correct on the Chevalley 
generators. For example, 

k 

The second line vanishes because the antipode is used to define the dual. Finally, (5) 
tells us the scalar factor for the inverse. Setting z — gives 

This completes the construction of change of basis matrices. 

5.5 Spontaneous polarisation. Two decades ago, Baxter calculated the spontaneous 
polarisation of the six- vertex model using the Bethe Ansatz [44]. Unlike many other 
important results, no alternative derivation was subsequently forthcoming, although 
the method has recently been extended to the calculation of the mean polarisation 
in the higher spin case [45]. The properties of VOs gives the first new derivation: 
moreover VOs give a general method to obtain correlation functions. The ideas have 
even been used for the eight- vertex model [46], validating a long-standing conjecture 
[47], although the underlying algebraic structure is not presently understood in that 
case. 

Here we sketch the derivation for the six- vertex model, which is rank-1. The 
solution for arbitrary rank has recently been given by Koyama [33] . The expectation 
value that a single edge variable has the value m (= ±) is simply the trace of an 
operator with matrix elements Ejk = Sjm^km- Recalling our identification of the 
CTM generator Kctm with p, consider the more general unnormalised expectation 
value 

F^\z,/Z2) = trv^(A^)(Q-2^$j!r;(^i)<^-^''(^2)). 

It takes its values inV* ®V and the spontaneous polarisation is the ratio (Fj*|(l) — 
F^l{l))/{F^l{l) + F^l{l)). 
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We may derive a g-difference equation for F*^*^: 

=tr^(^,)(g-^-$i;r:(-i)<^-^''(-2g^)) 

= -qz-'^^r{zq-'){id^q-^^)R%z)xtry^^^_^^{q-'P^^^^^^ 

The three steps use, respectively, the cychc property of the trace, the homogeneity 
and intertwining property of the VOs, and the commutation relations. The equation 
thus obtained reads 

F«(zg-^) = -qz-^''{id®q-^')rizq-^)R*{z)F^^-'\z). 

These equations may be solved quite simply. Set G^^\C) = F^^\C^)±qCF^^\C^), 
which decouples them. Note also that only the 2 by 2 block spanned by ® v± is 
non zero. For it we find 

This may be reduced to scalar equations by considering a generalised eigenvalue prob- 
lem: M{Qw{Q = X{Qw{(^q~'^). The eigenvectors are 

i^ = (i,±r'). A(o = c. 

The solution for G^^\Q involves 0^2 (±C), where 

0p(C) = (p;p)oo(C;p)oo(C~V;p)oo- 

It satisfies Qp{Cp) — — C~^©p(C)- order to have the required analyticity a zero of 
this function must cancel the factor (1=fC) '^±{0 c^nd this implies that there is only 
one non-zero component for each of G^^\Q. The argument whereby one shows that 
F^'^\z) has the assumed analyticity is based on finding an integral representation, 
and is given in [33,46]. The solutions may be normalised by the fact that F^'^\q~'^) 
is the trace of the identity operator. However, this is not necessary for finding the 
staggered polarisation. The G^^\Q differ only in the factors 0^2 (±C), and we find, 
independent of i, 



P = 



(i + C)e,2(c) 
(i-C)e,2(-c) 



Q (1 - q'"-^'? 

C=l n=0 



11 (1 _|_ ^2n+2-)2 
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The beautiful interplay which has emerged between physics and mathematics 
brings to mind a quote from the writings of Ludwig Boltzmann: 

It is unbelievable how simple and straight forward each result appears once 
it has been found, and how difficult it seems so long as the way which leads 
to it is unknown. 
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Appendix: Some numerical results 

For making finite size numerical calculations it is convenient to use -K'ctm rather 
than -f^xxz because the gauge term forces the boundary conditions to those appro- 
priate for V{Ao) when N is even and V{Ai) when N is odd [39]. Here we show some 
results for chains of length up to 11. Energies are relative to the antiferromagnetic 
ground state energy. They should be compared with the figures for the weight space 
structure of the relevant modules, t/g (312) -in variance means that the fi multiplets 
are exact. The fo multiplets tend to the proper degeneracy as N increases. Numerical 
checks have also been made on the actions of eo and fo for these values of N. In the 
weight space diagrams for the level-1 modules, spin is in the horizontal direction and 
grading is downwards. The action of fi is shown as a left pointing arrow, since /i 
decreases the weight by cti. /o is a diagonal arrow down and to the right, since fo 
decreases the weight by ao, and ao + ai = S. 



hwv 



Spill 


= 1 


= 6 


N = 8 


A^ = 10 

















1 


1.0021201 


1.0002473 


1.0000301 


1.0000038 





2.1189021 


2.0203449 


2.0033092 


2.0005036 


1 


2.1167936 


2.0200996 


2.0032794 


2.0004998 







3.2770532 


3.0678881 


3.0151265 


1 


3.1200455 


3.0205449 


3.0033418 


3.0005085 


1 




3.2768064 


3.0678580 


3.0151227 







4.2968794 


4.0711308 


4.0156211 









4.4835212 


4.1520062 


1 




4.2771017 


4.2771017 


4.0151267 


1 






4.4839116 


4.1520024 


2 




4.0205469 


4.0033421 


4.0005086 



Weight spaces of V{Ao) and negative eigenvalues of truncated CTM 
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hwv 



Spin 


iV = 5 


N = 7 


N = 9 


iV = 11 


1/2 














1/2 


1.0081308 


1.0011185 


1.0001513 


1.0000204 


1/2 


2.1899754 


2.0396425 


2.0076092 


2.0013141 


3/2 


2.0082121 


2.0011297 


2.0001528 


2.0000206 


1/2 


3.1980888 


3.0407572 


3.0077598 


3.0013343 


1/2 




3.3752382 


3.1052883 


3.0268467 


3/2 


3.1899870 


3.0396431 


3.0076093 


3.0013141 


1/2 


4.2000343 


4.0411612 


4.0078363 


4.0013475 


1/2 




4.3763556 


4.1054395 


4.0268671 


1/2 






4.6000763 


4.2077221 


3/2 


4.1999428 


4.0411497 


4.0078347 


4.0013473 


3/2 




4.3752396 


4.1052884 


4.0268467 



Weight spaces of V^(Ai) and negative eigenvalues of truncated CTM 
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